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In this paper we consider a generalization to the asexual version of the Penna model for biological 
aging, where we take a continuous time limit. The genotype associated to each individual is an inter- 
val of real numbers over which Dirac <5-functions are defined, representing genetically programmed 
diseases to be switched on at defined ages of the individual life. We discuss two different continuous 
limits for the evolution equation and two different mutation protocols, to be implemented during 
reproduction. Exact stationary solutions are obtained and scaling properties are discussed. 
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I. INTRODUCTION 

A possible explanation for aging in biological populations is based on deleterious mutation accumulation whose 



> 

• effects are felt at late ages when the intensity of natural selection is lower: mutations in the genome causing diseases 
genetically programmed to happen at late ages will not prevent reproduction and these mutations may spread through 
> the population [^]J^] . Population dynamics always involve non- linear evolution equations and the attractors of these 
dynamical systems may develop correlations in the stationary solutions that cannot be explained by the analysis of 
the equations per se. On the other hand, analytical solutions to these equations may be very hard to obtain. 

The Penna model for biological aging Q has been proposed as a tool to perform numerical simulations of dynamical 
systems describing age-structured populations subject to genetic mutations and to natural selection, allowing quanti- 
tative analysis of the possible attractors. It has been successfully applied to many different aspects of biological aging, 
considering asexual and sexual reproduction, haploid and diploid populations, different reproduction strategies, etc.. 
For recent reviews see [f|||]. 

The asexual version of the Penna model considers a bit-string associated to each individual of a given population. 
This bit-string is inherited from the individual's parent and contains information about the genetically programmed 
diseases that the individual will develop during its life. A one in the k th bit of the string implies a disease to be 
switched on at age k. Each individual dies whenever the T th disease happens in its life, that is, the locus of the T th 
bit set to one gives information about the inherited maximum age this individual may reach. To avoid an exponential 



growth of the population, a Verhulst factor is considered through a survival probability given by (1 — N(t)/N max ), 
where N(t) is the total population at time t and N max is the carrying capacity of the environment. At each time 
step the living individuals grow older one time unit - usually called year - and the corresponding bit in its bit-string 
is read, testing for the genetically programmed death. After a juvenile period, ending at age R, the individuals start 
to reproduce. Different assumptions on the fecundity after R have been considered to describe different situations 
ranging from constant fecundity up to age i?/, with no reproduction after that, to increasing with age fecundity, 
describing trees or crustaceans || that can grow very old. Also, different protocols have been proposed to produce 
the offspring bit-strings. Asexual reproduction is simulated by considering replicas of the parent bit-string except 
for M bits that may suffer deleterious mutations. Sexual reproduction protocols mix the information contained in 
both parents bit-strings and then consider mutations. There are many different results, for different parameters and 
protocols, but it is clear that in general the model grasps the central issue of decreasing natural selection intensity 
with age, yielding good results for the age structure of populations. Analytical approaches to these models have been 
considered for the asexual Q and sexual || versions of the Penna model, confirming simulations results. 

Many parameters must be considered in the models. Some of them are inescapable ones: there are many different 
systems and species in nature and if one wishes to describe each of them as a consequence of such a general law 
as natural selection, the model to be proposed should allow for different sets of parameters, as mutation rate, for 
example. However, some of the parameters are intrinsic to the numerical techniques used, such as the length B of the 
bit-strings. These parameters may be turned inoffensive if scaling could be found, that is, they could be related to 
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natural units of the system. In a recent paper, Malarz |9[] considered simulations of the Penna model with different 
bit-strings lengths and concluded that a scaling may not be possible. 

When scaling is to be considered, it is always interesting to avoid discrete equations, since the size of the time steps 
relatively to some parameters may be relevant to the scaling. Here we consider the asexual version of the Penna model, 
as proposed in reference Q, and take its continuous time limit. It turns out that age must also be continuous, and 
bit-strings turn into strings of real numbers of semi-infinite length, over which a sum of Dirac 5-functions is defined. 
We show that the lack of scaling is a consequence of the mutation protocol used, when the mutation probability is 
proportional to the inverse of the bit-string length, B . On the othar hand, when a Poisson distributed mutation 
probability is considered, as suggested by Bernardes and Stauffer [Q, we show that the results are independent of the 
string length considered, provided it is longer than the maximum possible age predicted for the population. We do 
not contradict previous work without scaling since that work |j] used a mutation probability that is proportional to 
the inverse of the bit-string length. Moreover, besides the scaling with the length of the bit-string, we also show that 
further scaling properties arise when an adequate continuous time limit is taken to the evolution equations, allowing 
the definition of R as a natural time unit for each population. 

The paper is organized as follows. In the next section we take the continuum limit of the discrete model for 
asexual populations by considering two possible continuous time evolutions. In section 3 we analyze possible mutation 
operators and in section 4 we present solutions and discuss scaling properties. Finally in section 5 we discuss the 
results and conclude that a full scaling may be found only when adequate mutation protocol and time evolution 
equations are considered. 



II. FROM DISCRETE TO CONTINUOUS TIME 



An analytical approach for the asexual, discrete version of the Penna model considers relative populations 
x(a,m,t), given as 

x(a, to, t) = N(a, m, t)/N max , (1) 

where N max is the carrying capacity of the environment and N(a,m,t) is the number of individuals at time t, with 
age a and death age to, that is, the age at which the T th bit one of its genome will be read. The time evolution of 
these populations is described by the following set of discrete time equations: 

x(a + 1, m, t + 1) = (1 — x(t)) x(a,m,t) for < a < to, . , 

x(a + 1, m, t + 1) = for a > m. 

Here x(t) is the total population at time t, that is, it is the sum over ages a and death ages m of x(a, m, t). The 
equation for a = states how offspring is produced. The related equation is 

x(0,m,t + l) = (l-x(t)) ^ ^2 F(m,m') x(a,m',t) , (3) 

a— R m'— a+1 

where the birth matrix F(m, to') is the probability that a parent with death age to' gives birth to a child with 
death age m. When only bad mutations are considered, F(m,m') is triangular and it is possible to obtain the exact 
stationary solution for the above evolution equations (tJ, provided the birth matrix fulfills the following conditions: 

1. F(m,m') is a triangular matrix such that F(m,m') — for to > to', that is, parents cannot give birth to 
offspring with larger life expectancy (it corresponds to only bad mutation in the Penna model). This is not 
biologically unrealistic for well adapted populations - advantageous mutations are expected to be extremely rare 
due to the large times required for noticeable species evolution; 

2. F(m, to) ^ 0, that is, the probability that the parent gives birth to offspring with its expected life length is 
different from zero. This condition is also expected in biological populations, and 

3. F(m, to) < F(m' ', to') if to > to', that is, the probability that a parent gives birth to offspring with its expected 
life length to decreases with to. In other words, the larger the parent expected life length, the larger the 
probability that a difference in the genetic charge of the offspring effectively reduces their expected life length 
(the better the genetic code, the larger the number of events that can spoil it). 
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In the original Penna model, mutations are implemented in the following way: for each newly produced bit-string, 
M loci on the bit-string with B bits are randomly chosen and set equal to 1, regardless their previous state. The 
birth matrix for the asexual Penna model can then be explicitly obtained by considering the probability of one 
mutation happening in a locus with bit zero, and estimating the change in death age caused by that mutation. The 
birth matrix thus obtained fulfills the above conditions. 

The results depend on the reproduction strategy. For a given initial reproduction age R and same mutation rate 
M, it can be proven that for final reproduction age Rf — > oo, with constant fecundity after R, there is a maximum 
possible life span m max for the stationary state. For initial populations containing m > m max , if R < m max the 
population evolves towards an age structure that gives a survival rate fairly constant up to R, and then decaying 
smoothly to zero asa-> m max . If, however, Rf < m max or the initial population has a maximum life span to* smaller 
than m max , the survival rate for the stationary solution decreases smoothly up to Rf or to*, where it abruptly decays 
to zero. In the special cases where Rf = R or R > m max , reproduction happens for only one year and the whole 
population die after it: this is the case of semelparous species that present catastrophic senescence [0J|. 

Here we present a continuous version of this model, and discuss different ways of implementing mutations. We 
start by considering time as a continuous variable, hence age should also be continuous. As death can happen at any 
moment in time, the genome must be described by a real variable interval. To each individual we associate a genome 
function f(£) of a real variable < t < L, where the string length L may possibly go to infinity. This function is a 
sum of Dirac (5-functions, that is 

= -4), (4) 

»=i 

where the locations of the (5-functions, £i, represent the ages at which genetical diseases are switched on. We define 
the health status of an individual as 

H(a)= f f{l)dt (5) 
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such that the death age of the individual is reached when H(a) — T, that is, when the T th (5-function - or the T th 
genetically programmed disease - happens in its life. Now x(a,m,t) da dm is the relative number of individuals at 
time t with age between a and a + da and death age between to and m + dm. The total population x{t) is calculated 
from x(a, m, t) as 



x(t) = / dm / da x(a,m,t) , (6) 



where m* is the maximum death age in the population (to* < m max ). There are two continuous time evolution 
equations for x(a, m, t) that are compatible with the discrete limit Eqs. (0), as we will discuss now. 



A. The Volterra description 

This approach is based on the product integration calculus proposed by Volterra in 1887 [ pT| , ^2[ . The first of Eqs.(|^) 
can be written in the continuous version as 

x(a + da,m,t + dt) \ 1/{Kdt) 

xja^rnj) ) = (1 " * (<)) ' (7) 

where the left hand side has been taken to the power 1/ (ndt) and k is a constant with dimensions of time -1 such that 
nAt = 1 in the discrete limit: k _1 is a convenient time unit. The above equation can be rewritten as 

— lna;(a, m, t) = /tln(l — x(t)) , (8) 

where the total derivative can be displayed as 

d , , d\nx{a,m,t) da d\\ix(a,m,i) dm d\nx[a,m,t\ 
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Clearly da/dt = 1 and dm/dt = 0. The stationary solution ( dln x ^ m ' t ) = o) for Eq.(||) is then 

x(a,m) = x(0,m)(l-x) Ka , (10) 

where we have dropped the time index for the stationary solutions. Observe that this is in the same form as the 
solution for the discrete version of the model ItI . 



B. The orthodox description 

A more orthodox way of writing the continuous limit for Eqs. (|2|) is to rewrite them as follows 

x(a + da, m, t + dt) — x(a, m, t) 

= - x {t) x{a, m, t) , (11) 

ndt 

where k has the same role as in the Volterra description: a convenient unit such that for the discrete limit time 
interval nAt = 1. The differential equation is 

—x(a,m,t) = — k x(t) x(a,m,t) , (12) 
dt 

with the total derivative analogously given as in eq. (|9|). The stationary solution is 

x(a,m,t) = x(0, m) exp (— nxa) . (13) 

The difference between stationary solutions Eqs.(|l0|) and ( |T3] ) is not important when x ~ 0, since in this case 
(1 — x) Ka ~ exp (— Kxa). On the other hand, observe that the above equations allow x > 1. Both solutions depend on 
the newborn populations x(0,m) that in turn depend on the birth matrix F(m,m'), that is 

pm* rrn 

x(0,m)=b / dm' / F(m,m') x(a,m') da, (14) 
Jr. Jr. 

where b is the number of offspring produced per fertile individual per time unit. 
Before proceeding with the solution let us obtain F(m,m'). 



III. THE BIRTH MATRIX AND MUTATION OPERATORS 



Here we shall consider two ways of implementing mutations, an analogue to the discrete Penna model and a Poisson 
distributed probability of happening mutations. 



A. Penna model mutations 



Let us assume that the genome length is L. One parameter of the model is the number M of loci randomly chosen 
to assign a new 5-function on the genome. Diversely from the discrete model, the probability of this happening on 
a locus already with a 5-function is vanishingly small. Hence, the only possibility of the offspring death age being 
equal to that of the parent occurs when no new mutation happens before mf. The mutation operator A{m,m!) for 
exactly one mutation is given as 

Aim, m') = h—UL S( m - m >) + Z( — f-i gr m > _ m \ r 15 \ 
L L m 

where 6(m! — m) is the step function (equals to 1 if m! > m and otherwise) that guarantees that the child death 
age is always less or equal to that of the parent. The birth matrix for M mutation is 

F(m,m') = (A(m,m')) M . (16) 

^From Eqs.(^|) and ( |l6| ) it is clear the origin of the lack of scaling: as L is varied, the mutation probability per unit 
of genome length varies. And as T appears both as a linear coefficient and as an exponent, it is at least not trivial 
to find sets of parameters that would yield equivalent solutions. On the other hand this mutation protocol is too 
artificial: certainly there are genomes that do not suffer mutations at all as there could be some that suffer more than 
a fixed number of mutations. A more realistic assumption is to consider a fixed (small) mutation probability per unit 
of genome length, taken from a Poisson distribution, as we consider in what follows. 
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B. Poisson distributed mutations 



Consider the probability of n mutations happening in a given length i of a genome as given by a Poisson distribution 

(pq-expHrf) 
n\ 

where p is the probability of occurring one mutation (adding a new 5-function) per unit of genome length. To build 
up the birth matrix F(m,m'), there are three possibilities to consider: i) if m > ml ', then F(m,m') = 0, since we 
consider only the possibility of adding new 8— functions to the genome (deleterious mutations); ii) if m = m', meaning 
that no mutation occurs before ml and the offspring death age is the same as that of its parent, and finally Hi) if 
m < to', when new 6— functions are added to the genome before the occurrence of the parent T th disease. Taking all 
the possibilities into account, the birth matrix can be written as 



F(m, m) — S(m — m) exp (—pm') 



to) 



to k to -to T ^_ fe {pm) 1 



to' to' (T — 1 — k)\ 



T-2 , , \T-\-k 



+ cxp( _ pm) y^. 2 (^) fc { m '~ m ^-2- k (pm) 



'to' to' (T — 1 — k)\ 



(18) 



where 6 {ml — to) guarantees that F(m, to') — if to > to' and C T = T!/(fc! (T - fc)!). Observe that mutational 
meltdown may be prevented since there is a non-zero probability that some offspring are bred without additional 
harmful mutations (to = to') p0[ . It can be shown that the above birth matrix obeys the normalization condition 

dm F{m,m) — 1 , (19) 



o 

that is the sum over all possibilities for the offspring genome. Comparing Eqs. ( |l5| ) and (|lj), we observe that the 
mutation controlling parameters in the different mutation protocols are respectively L and p, that is, using Poisson 
protocol the mutation probability decouples from the genome length L. Hence, since L no longer appears in the 
mutation operator, the results do not depend on L, and scaling should be expected. When only bad mutations are 
considered, there may exist a maximum life span m max to the steady state solutions that is independent of L. In 
this case, the stationary solutions are not affected by bits located after m max and consequently are independent of 
L, provided the mutation rate per genome length does not depend on L. In other words, the dependence on L seems 
to be responsible for the lack of scaling in Eq.(|l5|) when mutations are implemented as in the discrete version of the 
Penna model. 

On the other hand, further scaling properties can be found in the solutions of the evolution equations, provided the 
Poisson mutation protocol is used, as we show in the next section. 



IV. SOLUTIONS AND SCALING 



To have further insight of what happens we shall consider explicit solutions to the evolution equations. From now 
on we shall consider only Poisson mutation protocol. We first rewrite Eq.(p^) for the newborn offspring using Eq.(|l3|) 
for the orthodox evolution: 

b f 00 

x(0,m) = — / dm' [exp (— xkR) — exp (— XKm)] F(m,m') x(0, m) . (20) 

We observe next that x(a,m) has dimensions of time -2 , 6, p, and F(m,m') have dimensions of time -1 , while to, R, 
and a have dimensions of time 1 and x is dimensionless. We can take k~ 1 as the time unit and rewrite all equations 
in dimensionless variables, which is equivalent to rewrite all equations setting k = 1 and reading all quantities as 
given in time unities of k . We then have one less parameter, and we remark that reading p in units of k links the 
mutation rate p to the exponential constant of time evolution. 

Even with one less parameter, Eq.(|2^) is an integral equation, involving all newborn population and the birth 
matrix. We can solve this equation because F(m, to') is triangular, but we must find out the maximum possible 
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value m max for the death age that any stationary solution may present. We do this in an analogous fashion to that 
used for the discrete version of the model R] . We assume that there is a maximum death age to* in the population. 
As F(m,m') is triangular - there are only bad mutations - as time proceeds this maximum value cannot increase, 
although too long living individuals may disappear. Offspring with death age equal to m* can be produced only by 
m*-parents. Using Eq. (|l8|) , the newborn equation for to = to* can be written as 

x(0, to*) = — [exp (— xR) — exp (— xm*)] exp(— pm*) x(0,m*) . (21) 

The trivial solution, x(0, to*) = 0, is always possible. A nontrivial solution, x(0, to*) =^ 0, is possible if two conditions 
are met. First, 

— exp (— pm*)[exp (—xR) — exp (—xm*)] = 1 , (22) 
x 

what guarantees that Eq.(^i"|) is satisfied, and a second condition that guarantees that x(0,m* — dm) > 0, 

— — exp (— pm)[exp (—xR) — exp (—2:771)] | m=m » > , (23) 



777* < 77w = R + - In ( V -^~ ) • (24) 



resulting in 

777* < m max = R + 

x \ P 

All these calculations may be repeated for the Volterra evolution equations. The equivalent equations to be solved 
have the same form as Eqs.(p2[) and ( p3|) with x replaced by — In (1 — x), and the corresponding limit for the maximum 
life span in a given population is 

m* < m max =R - - 1 - In ( P - Hl - X) ) . (25) 
ln(l - x) V P ) 

A limit for the life span of a population means that if, at initial times, the maximum life span in the population is 
777,* > m max , in the stationary solution the population goes to a state where m* = m max . But if initially to* < m maxi 
the population remains with that maximum life span, since in this model only bad mutations are allowed to happen. 
For each possible to* (< m max ) the total population x may be obtained from Eq.(p2[) or from its analogue for 
respectively orthodox or Volterra evolutions. m max is the largest value of to* that satisfies both conditions ( p2| ) and 
( p3| ) or their analogues. 

These two limits for the maximum life span for different evolution protocols depend differently on the value of the 
total population. We can observe some scaling properties regarding the maximum possible life span m max and total 
population x. If we have, for every a > 0, 

p' = ap , 

R' = -, (26) 
a 



Then for the orthodox evolution 



b' = a b . 



x = a x , 

7 = 

max 

a 

„3„ 



(27) 



x'(a/a,m/a) = a 3 x(a,m,) . 
On the other hand, for the Volterra description the total population scaling is different: 



x' = 1 - (1 - x) a , 

(28) 



a 

and the relation between population densities a;(a,7T7) and x f (a/a,m/a) is far from trivial: 



G 



1— / dm / da x(a,rri) = 1—1 d(m/a) / d(a/a) x'(a/a,m/a) . (29) 
■/o Jo J Jq Jo 

We note that for both evolutions the value of maximum life span does not depend on T. Hence varying the number 
of harmful diseases that the individuals may tolerate will not change the maximum life span, although it can change 
the survival rate for a < m max - The reason is that both conditions stated in Eqs. ( p2| ) and ( p3| ) are generated by the 
first term in the expression for the birth matrix, Eq.(p"8|), that gives the probability for a child with the same death 
age as its parent. This term depends only on the parent death age, and not on the number of diseases that there are 
in the parent genome before to'. This is also in agreement with the fact that the increasingly better medical care has 
increased human life expectancy but has not equally changed the maximum life span |l3| , |l4| ] . In the Penna model the 
effect of proper medical care is equivalent to a shift of the maximum number T of genetically programmed diseases 
that an individual may endure p5| . In Fig.(|l|) the plots of maximum life span and total population as functions of 
mutation probability p for different values of fecundity b. The maximum life span decreases with increasing p and 
decreasing b. 



£ 
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FIG. 1. Plots of (a) maximum life span and (b) total population as functions of fecundity b for different values of mutation 
probability p. The maximum life span decreases with increasing p and decreasing b. Observe that here x can be greater than 
one, since we are considering the orthodox time evolution. 

To fully appreciate the scaling present in the solutions we must solve Eq. (p0|) . For that we first transform the 
integral equation into a differential one. Consider a stationary population with to* = m max , then for m > R the birth 
equation can be rewritten as 

exp (-xR) - exp (-xm) \ 

l — o exp (— pm) x(0, m) = 

x J 

u i \i m \T-\ [ mmax ex P {-xR) - exp {-xm') , . . , 

b p exp (— pm)( — -J / j[m - m) x{0,m) dm , (30) 

771 J rn X 

where /(to' — to) is defined as 



f(m' ~ m) = E C^^-^ + J2 C$ty (3D 



k=0 k=0 



such that, after some calculations we can write 

F(m, m!) = exp (—pm) 5{w! — to) + p exp (~pm) ( — -) T ^ 1 f(m' — m) 9(m' — m) . (32) 



We now rewrite Eq.(30) as 



exp (-pro) A . - exp(-xR) -exp(-xm)\ 

'1 — 6 exp {—pm) 1 x[0, to) 



b m T 1 V i- \ jt / ^ 

" exp (— xR) — exp (—xm') 



f(m — ro) x(0, m) dm' . (33) 
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We then differentiate T times both sides, and use the fact that we can caiculate the value of f(rn' 
derivatives at m = m! . We then arrive to the following differential equation: 



d T ( exp (pm) 



x(0, m) I = exp (pm) 



d T ( exp (—xR) — exp (- 



dm T 



,T-1 



exp (—pm) x(0, m) 



m) and all its 



(34) 



This equation may be numerically solved. To obtain the analogue of the above equation for the Volterra evolution it 
suffices to replace x by — In (1 — x). Nevertheless, in what follows we will consider only the orthodox evolution. 
We have numerically solved Eqs.(|34]) for T = 1 and 
T = 2. Due to the scaling properties we can use the ini- 
tial reproduction age R as our 'natural' unit. In Fig. (2) 
where we plot x(0, m)/x(0, R) for sets of R, p and b re- 
specting the scaling conditions given in Eqs.(|26]), for 
both T = 1 and 2. The scaling is clearly shown. In 
what follows we shall use R = 1 and mutation and 
birth rates p and b are given in units of RT 1 . 

The effect of different mutation and birth rates on 
the population can be monitored through the survival 
rate S(a) that is defined as 



S(a) = 



x(a + da) 



l/da 



exp( 



x(a) 

dhi(x(a)) , 
da 



(35) 



Figure (3) shows how the survival rate changes when 
the fecundity b varies for the same value of mutation 
rate p — 0.001. As b increases, the total population x 
also increases and the survival probability is uniformly 
reduced at all ages. Also, the maximum life span and 
the onset of senescence, marked by the strong reduc- 
tion in survival rates, are shifted to earlier ages. The 
figure shows data for both T = 1 and 2. The maximum 
life span m max remains the same for different T, but 
one can see that at the same age, higher values of T 
imply lower survival rates, that is, better medical care 
decreases the strength of natural selection such that 
more mutations accumulate at later ages. 

Figure (4) shows how the survival rate depends 
on the mutation rate p, for the same value of b = 4. 
As expected, a smaller p enhances the maximum life 
span. However, young individuals present lower sur- 
vival rates due to larger total population x, that re- 
duces mortality through the exponential pre-factor 
cxp(— xa). 

Finally, Fig. (5) presents the normalized mortality 
function q(a), defined as 



q(a) = — In 



S(a) 
S(0) 



(36) 



for different values of 6, p and T. A Gompertz region is 
clearly seen, and also a deviation from the exponential 
behaviour for young ages. 




FIG. 2. Plots of x(0,m)/x{0,R) for sets of R, p and 6 re- 
specting the scaling conditions, as given in the inset, for both 
T = 1 and 2. This behaviour refers to the orthodox evolution 
with Poisson mutation protocol. 




FIG. 3. Survival rate S(a) versus b for mutation rate 
p — 0.001 for the orthodox evolution with Poisson muta- 
tion protocol. Observe that smaller T imply higher S(a) for 
the same a. 




FIG. 4. Survival rate S(a) versus p for birth rate b — 4 for 
the orthodox evolution with Poisson mutation protocol. 
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V. CONCLUSION AND DISCUSSION 



In this paper we have generalized to continuous 
time the discrete asexual Penna model for aging, such 
that the genome of the individuals is now represented 
by a sum of 8- functions on continuous interval, 8 (a — 
do), each one representing a disease to be switched on 
at the age a^. 

Two different time evolution protocols could be 
used: Volterra and orthodox time evolution equa- 
tions are compatible with the discrete formulation of 
the Penna model. However, we showed that a com- 
plete scaling is only found with the orthodox evolution 
equations (eqs. @-@ ). 

Also, two different protocols may be considered to 
describe mutations happening in offspring genomes. 
One considers a mutation rate that is proportional to 
the length of the genome. This is equivalent to the 
mutation protocol commonly used with the discrete 
version of the Penna model, where scaling may not be found due to the somewhat artificial dependence of the 
mutation rate on the length of the genome. A second mutation protocol considers a Poisson probability distribution, 
with mutations occurring with a fixed probability p per unit of genome length. As this last mutation protocol does 
not depend on the genome length, the whole set of stationary solutions does not depend on L. 

Morever, in the special case of the orthodox time evolution and Poisson mutation protocol, we could find further 
scaling relations such that the initial reproduction age R may be taken as the natural time unit of the system. In this 
sense, all rates, as mutation and birth rates, are renormalized in units of R . One direct consequence of this scaling 
is the correlation between earlier initial reproduction age and earlier senescence, an effect that is due to the scaling 
inherent to the model that shows up in the stationary solutions of the population dynamics equations. That is, here 
demographic effects only are responsible for this correlation and antagonistic pleiotropy effects are not needed. 

The generalization of sexual version of the Penna model is now under investigation and shall appear elsewhere. 
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